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(54) Magnetic resonance imaging apparatus 

(57) A magnetic resonance imaging apparatus in- 
cludes superconducting magnets (10) which generate a 
temporally constant magnetic field through a bore (12). 
The bore, hence the superconducting magnets, have a 
length which is relatively short compared to its diameter, 
eg a length to diameter ratio of less than 1 .5 to 1 .0 and 
typically 1 :1 . A gradient coil assembly (30) of the appa- 
ratus is disposed around the bore for generating gradient 
magnetic fields across the bore. With such a relatively 
short bore magnet, in the region in which the gradient 
field coil is disposed, the main magnetic field suffers 
non-uniformities including radial magnetic field compo- 
nents. When current pulses are fed through the windings 
of a primary gradient coil (32) and a secondary gradient 
coil (34) of the gradient coil assembly the currents inter- 
act in an unbalanced manner with the non-uniformities 
and radial components of the temporally constant mag- 
netic field, causing a net force in axial and/or transverse 
directions. Force correction windings (64, 72) are mount- 
ed adjacent the ends of the gradient coil, carrying current 
in an opposite direction to z and transverse gradient 
windings in order to produce a counterbalancing force 
such that the net and counterbalancing forces substan- 
tially cancel. In this manner, the gradient coil assembly 
is freed from the necessity of a substantial mechanical 
mountings and the vibration and distortion which large 
net forces can cause. 
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Description 

This invention relates to magnetic resonance imaging apparatus. It relates particularly to magnetic resonance im- 
aging apparatus employing "short" bore magnets and will be described with particular reference thereto. However, the 

s present invention finds application in any magnetic imaging apparatus wherein a gradient coil is disposed in a magnetic 
field such that Lorentzian force components on the coil are not balanced. 

Magnetic resonance imagers commonly include a bore having a diameter of 90 centimetres or more for receiving 
the body of an examined patient. The bore is surrounded by a series of annular superconducting magnets for creating 
a substantially uniform magnetic field longitudinally along the patient receiving bore. The more axially spaced the annular 

10 magnets, the more uniform the primary magnetic field within the patient receiving bore tends to be and the longer the 
axial dimension over which such magnetic field uniformity exists. Typically, the bore is at least 1 .6 metres long and often 
longer. 

One of the drawbacks to such "long" bore magnets is that the region of interest is often inaccessible to medical 
personnel. If a procedure is to be performed based on the image, the patient must be removed from the bore before the 
15 procedure can be performed. Moving the patient risks potential misregistration problems between the image and the 
patient. 

One way to improve access to the patient is to shorten the length of the magnet and the patient receiving bore. If 
the magnet and the bore were shortened to about 1 metre or roughly the diameter of the bore, patient access is much 
improved. Although the size of the uniform magnetic field area compresses to a more disk-like shape, the area of sub- 

20 stantial uniformity is still sufficient for a series of 1 0 to 20 contiguous slice images. NMR helical or continuous scanning 
methods can also be employed. 

Although an adequate imaging volume remains, the magnetic field in the volume around the periphery of the bore 
which receives the gradient coil tends to become relatively non-uniform and has both axial z-components and radial 
x,y-components. The gradient coil generally includes windings for generating three linear and orthogonal magnetic field 

25 gradients for providing spatial resolution and discrimination of nuclear magnetic resonance signals. Gradient coils are 
typically designed and constructed to optimize strength and linearity over the imaging volume and stored energy and 
inductance in the gradient coil. See, for example, U.S. Patent No. 5,296,810 of Morich. To create the magnetic field 
gradients, current pulses are applied to the x, y, and/or z-gradient coils. These currents interact with the main magnetic 
field to generate Lorentz forces on the gradient coil. Due to the symmetries included in gradient coil currents, the Lorentz 

30 forces across the entire coil cancel when the gradient coil is disposed in a uniform magnetic field. However, when the 
main magnetic field is less uniform, particularly when there are significant radial and non-uniform axial components in 
the neighborhood of the gradient coils, a net thrust can be developed. Typically, pulsing the z-gradient coil causes a net 
thrust in the z-direction, pulsing the x-gradient coil develops a net thrust in the x-direction, and pulsing the y-gradient 
coil causes a net thrust in the y-direction. In the case of the z-gradient coil, the net axial force can be on the order of a 

35 few hundred pounds. These net thrusts tend to push or urge the gradient coil axially out of the bore. Although the gradient 
coil can be anchored mechanically, these large forces still tend to cause acoustic noise and increased vibrations to the 
magnet and the patient. Such vibration has deleterious effects on imaging, such as a loss of resolution. 

The present invention provides a magnetic resonance imaging apparatus including a gradient coil which alleviates 
the above-referenced problems. 

40 According to the invention there is provided a magnetic resonance imaging apparatus which includes a main magnet 

for generating a temporally constant magnetic field through a patient receiving bore and a gradient field coil assembly 
and an RF coil around the bore and in which the temporally constant magnetic field has components which interact with 
current pulses applied to principal windings of the gradient field coil assembly to generate-magnetic field gradients in 
such a manner that a net force is experienced by the gradient field coil assembly, characterised by force correction 

45 windings mounted adjacent edge areas of the gradient field coil assembly which carry currents in an opposite sense to 
the principal windings such that the net force is substantially cancelled. 

In an apparatus according to the invention the patient receiving bore typically has a length to a diameter ratio of 
less than 1.5:1, e.g. of substantially 1:1. 

In one particular embodiment of the invention the gradient coil assembly is a self-shielded gradient coil having a 

50 primary coil and a shield coil, current pulses passing through the primary and shield coils generating magnetic fields 
that combine within the gradient coil assembly to form a linear magnetic field gradient and tending to cancel outside the 
gradient coil assembly, and the current pulses that pass through the shield coil interact with radial magnetic field com- 
ponents to generate a first force, the primary gradient coil being wound such that the current pulses passing through 
the primary gradient coil interact with the radial magnetic field components to generate a second force that is substantially 

55 equal and opposite to the first force. 

In another particular embodiment the gradient coil assembly includes primary and secondary gradient coils which 
each include four symmetrically arranged thumbprint coils for generating a magnetic field gradient transverse to a lon- 
gitudinal axis of the bore, the thumbprint coils further generating a net transverse force, the force correction windings 
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including force offsetting current loops disposed adjacent an end of at least one of the primary and secondary gradient 
coils for carrying current flows which circulate in an opposite direction to a most adjacent thumbprint coil for offsetting 
the net transverse force along the transverse direction. 

In a further particular embodiment the gradient coil assembly includes a primary z-gradient coil including a series 
s of loop coils which are mounted around a cylindrical former for generating a z-gradient field and the force correction 
windings include loops extending around ends of the cylindrical former and connected in series with the z-gradient loop 
coils for carrying common current pulses but in an opposite direction. 

In a still further embodiment the gradient coil assembly includes a primary z-gradient coil and a secondary, shield 
z-gradient coil and the secondary z-gradient coil is disposed surrounding the primary z-gradient coil, current pulses 
10 applied to the primary and secondary z-gradient coils causing magnetic fields which combine within the gradient coil 
assembly to generate a linear gradient axially along the coil assembly and which generally cancel exterior to the coil 
assembly. 

The present invention also provides a method of designing gradient coils for a magnetic resonance imaging appa- 
ratus in which coil windings patterns are constrained to have a preselected gradient field strength, at least a preselected 
15 gradient field linearity and at least a preselected gradient field uniformity, characterised by: constraining the coil windings 
such that the gradient coil assembly experiences substantially no axial or transverse force when current pulses flowing 
through the gradient coil windings interact with radial and axial components of a main magnetic field within which the 
gradient coil is disposed. 

One advantage of the present invention is that it facilitates access to portions of the patient in the examination region. 
20 Another advantage of the present invention is that it improves image quality. 

Another advantage of the present invention is that it simplifies mounting and construction of gradient coils. 

One magnetic resonance apparatus in accordance with the invention and a method in accordance with the invention 
of designing the coils of a gradient coil assembly of the apparatus will now be described, by way of example, with 
reference to the accompanying drawings in which: 

25 

FIGURE 1 is a diagrammatic illustration of the magnetic resonance imaging apparatus. 

FIGURES 2A and 2B illustrate exemplary current distributions for z-primary and shield gradient coils of the appa- 
ratus; and, 

30 

FIGURES 3A and 3B illustrate a primary and shield or secondary x or y-gradient coil winding of the apparatus. 

Referring to FIGURE 1, the apparatus includes a plurality of primary magnet coils 10 which generate a temporally 
constant magnetic field along a longitudinal or z-axis of a central bore 12. The primary magnet coils are supported by 

35 a former 14 and received in a toroidal helium vessel or can 16. The can 16 is filled with liquid helium to maintain the 
primary magnet coils at superconducting temperatures. The can 16 is surrounded by one or more cold shields 18 which 
are supported in a vacuum dewar 20. 

A whole body gradient coil assembly 30 includes x, y, and z-coils mounted around the bore 12. Preferably, the 
gradient coil assembly is a self-shielded gradient coil assembly that includes primary x, y, and z-coil assemblies 32 

40 potted in a dielectric former 34 and a secondary or shield gradient coil assembly 36 that is supported on a bore defining 
cylinder 38 of the vacuum dewar 20. The dielectric former 34 with potted gradient coils can function as a bore liner or 
a cosmetic bore liner can be inserted to line it. Preferably, shims (not shown) for adjusting the magnetic field are positioned 
as needed between the primary and shield coil dielectric formers as taught by DeMeester, et al. in U.S. Patent No. 
5,349,297. A whole body RF coil 40 is mounted inside the gradient coil assembly 30. A whole body RF shield 42, e.g. 

45 a layer of copper mesh, is mounted between RF coil 40 and the gradient coil assembly 30. 

An operator interface and control station 50 includes a human-readable display such as a video monitor 52 and an 
operator input means including a keyboard 54 and a mouse 56. Track balls, light pens, and other operator input devices 
are also contemplated. Computer racks 58 hold a magnetic resonance sequence memory and controller, a reconstruction 
processor, and other computer hardware and software for controlling the radio frequency coil 40 and the gradient coil 

50 30 to implement any of a multiplicity of conventional magnetic resonance imaging sequences, including echo-planar, 
echo-volume, spin echo, and other imaging sequences. Echo-planar and echo-volume imaging sequences are charac- 
terized by short data acquisition times and high gradient strengths and slew rates. The racks 58 also hold a digital 
transmitter for providing RF excitation and resonance manipulation signals to the RF coil and a digital receiver for re- 
ceiving and demodulating magnetic resonance signals. An array processor and associated software reconstruct the 

55 received magnetic resonance signals into an image representation which is stored in computer memory, on disk, or in 
other recording media. A video processor selectively extracts portions of the stored reconstructed image representation 
and formats the data for display by the video monitor 52. An image printer can also be provided for making paper copies 
of selected images. 
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In the preferred embodiment, the diameter and length of the bore 12 have a size ratio of about 1:1. However, it is 
to be appreciated that the present invention is also applicable to other magnets, particularly magnets with limited main 
magnetic field uniformity. Typically, one might expect magnets with a bore length to diameter ratio of up to 1.5:1 to be 
candidates for the present invention. However, the present invention will also be applicable to magnetic resonance 

s imagers with longer bores in which there is a sufficient main magnetic field non-uniformity, particularly radial field com- 
ponents, in the area of the gradient coils to cause image degradation due to vibration or in which there are sufficient net 
forces that mechanical mounting of the gradient coils becomes difficult. 

The main magnetic field magnets 10 produce a main magnetic field B 0 over the imaging volume. In the short bore 
magnets, the z-component of the main magnetic field B z0 remains constant over a central imaging region of about 40 

10 to 45 centimetres in diameter in a x,y plane at the magnetic isocenter. Along the z-axis, the imaging region is shorter 
than 40 centimetres. Outside the ellipsoidal, uniform, imaging region, the B field changes significantly at larger axial or 
z-displacements from isocenter and with radial or p displacement from the central axis of the bore. The dependence of 
B z0 on the axial and radial positions becomes more significant around the gradient and RF coil locations, typically the 
region bounded by 30 cm. ^ p ^ 40 cm. and -60 cm. ^ z ^ 60 cm and encompassing all angular positions. Although 

15 inside the imaging region, the contribution of the radial component of the main magnetic field B p0 to the total value of 
the B 0 field is negligible, it is much more significant at the gradient coil sets' location, particularly in short bore magnet 
designs. 

The interaction between the spatially varying components of the magnetic fields and the current densities of the 
coil set creates two distinct problems. First, there is an interaction between the azimuthally directed current of the z-gra- 

20 dient coil with B p0 . Because both quantities are odd symmetric functions around the geometric center of the magnet 
plus gradient system, they generate a z-directed thrust force according to the Lorentz force equation. Depending on the 
axial or z-behavior of B p0 , in particular its value at the location of the conductors of the z-gradient coil, the value of the 
thrust force can reach several hundred pounds or a few thousand Newtons. 

The second problem deals with the interaction between the current density for a transverse x or y-gradient coil and 

25 both the B p0 and B z0 components of the main magnetic field. The result of such an interaction is again determined from 
the Lorentz force equation, but is a radially directed net force. Due to the azimuthal dependence of the transverse 
gradient's current density, the radially directed force is along the x or y-direction for an x or y-gradient coil, respectively. 
Again, the magnitude of this force is dependent on the value of the static magnetic fields B p0 and B z0 components at 
the locations of the current lines for the x and y-gradient coils. Because the radius of the x and y-gradient coils is different, 

30 the value of the corresponding force will also be different. The value of the components of B0 change with radial position. 

Looking first to the analytical evaluation of the z-directed thrust force which is the result of the interaction between 
the z-gradient coil and the B p0 component, a minimization technique is presented for designing z-gradient coils with a 
zero net thrust force. For a conventional z-gradient coil, the current is odd-symmetric around the geometric center of 
the coil and the magnetic field B p0 is also odd-symmetric around the geometric center of the magnet. Because these 

35 two centers normally coincide, the result of the interaction between the current and the B p0 field component is a z-directed 
thrust force on the gradient coil. 

With reference to FIGURES 2A and 2B, the z-gradient coil is self-shielded, having an inner or primary winding 60 
of radius a and an outer or secondary shield coil 62 of radius b. For simplicity of construction, the primary z-gradient coil 
is preferably circular loops wound on the dielectric former 34 and potted in epoxy or other resin. The shield z-gradient 

40 coil is preferably circularly wound on the dielectric former 38 that is incorporated into the vacuum dewar. Of course, 
alternative embodiments where it is not an integral part of the vacuum dewar are also possible. The current density for 
the inner gradient coil J^z) is azimuthally directed and varies only along the axial direction of the coil. Similarly, the 
current density for the outer gradient coil J (() b (z) is azimuthally directed and also varies along the axial direction of the 
coil. The interaction between these two current densities with the radial component of the main magnetic field at the 

45 radial locations of the coil, i.e., B p0 a and B p0 b , respectively, generates a Lorentz force which is directed along the z-di- 
rection. 

In general, the Lorentz equation describes the elemental force experienced by a current element J (() a ' b (z)dz(a,b) 
dcjja^ in the presence of the magnetic field B p0 a ' b as: 

50 df = D» (z)dzad<K^x B ^<J p + j^(z)dzbd«|x^xB^<jl p (1). 

The net force which is generated from these The net force which is generated from these two coils is along the 
z-direction and has the form: 
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with 



J, net 



4 



EP0 701 143 A1 



F z net = "2n ja f~j* (z) B p a 0 (z) dz +b f~jf (z) B p b 0 (z) dz | (2) , 
where B p0 a ' b (z) are known real functions of (z). The Fourier transform pairs for \^ a ' b and B p0 a ' b are: 



. a,b 
J* 


(k) 


= /_fj„ a ' b (z) e ikz dz 


(3a) , 


• a,b 


(z) 


= f~j*' b (k) e* ikz dk 
2n J-« 


(3b) , 


n a,b 


(k) 




(3c) , 


B P0 


(z) 


= ^/3 a o' b (k)e +ikz dk 


(3d) . 



These equations can be used to obtain the Fourier or k-space representation of the net thrust force as: 

F z net ={a f~lf (k) B p a 0 (k) dk + b f_+Jj* (z) B( £ (k) dkj (4) . 

Because B p0 a ' b (z) is a pure real, odd-symmetric function of z: 

B p o <a ' b> (k) = /~B p a 0 ' b (z)e^dz = -B p a 0 ' b (k) (5) . 

For a set of z-gradient coils with radius a,b, respectively, the expression of the axial z-component of the magnetic 
field due to the gradient coil B z , in the three regions which the two coils define can be simplified, due to the behavior of 
the current density. More specifically, no angular dependence is considered for a current which is uniform around the 
circumference of a cylindrical gradient coil. The current in the z-gradient coil is restricted to vary along the z-direction. 
With these constraints, the expression of B z in the three regions is: 

B z = 

^/"dke^kjaj^ (k) K, (ka) I 0 (kp) + bj 0 b (k) K, (kb) I Q (kp)] 
for p<a (6a) , 



^f-dk 
2nJ- 



e ikz k 



-a (k) K Q (kp) I x (ka) +bj # b (k) K a (kb) l 0 (kp) 
for a<p<b (6b) , 



bj^ (k) K 0 (kp) i x (kb) 
(6c) , 

where j ( j ) a (k), j ({) b (k) are the Fourier transforms of the current densities of the inner coil of radius a and the outer coil of 
radius b, respectively. Ideally, the z-gradient coil generated magnetic field B z is linear in the internal region of both coils 
and zero outside the two coils. In order to satisfy the constraint that the field is zero outside of the two coils, B z in the 
region b<p is set to zero, i.e., Equation (6c) is set to zero. One way to set Equation (6c) to zero is by relating the current 
densities of the inner and outer coils such that their Fourier components sum to zero, i.e.: 



B 



2n 



ikzi 



aj«f (k> K 0 (kp) i 1 (ka) - 
for b<p 
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or: 



<(k) K 0 (kp) I, (ka) + bj $ ° (k) K 0 (kp) I, (kb) = 0 
al, (ka) 



C (k) = 



(kb) J * 



j* ( k ) 



(7). 



(8). 



Substituting Equation (8) into Equation (6a), the expression for the magnetic field for the region inside of the two coils 
p<a becomes: 



B = 



2n 



I 0 (kp) Kj (ka) 



J^dkke^j^ (k) 
I 1 (ka)K 1 {kb) 



Z t (kb) K x (ka) 
From Equation (8), the net force expression can also be simplified: 



(9) 



= a 1 



/~j<f (k) | B* (k> --El^B* (k) 



dk (10). 



The stored energy in the coil W in terms of Equation (8) is: 



w=^£dk|j, a (k)r 



I 1 (ka) jq (ka) 



Ij (ka) K x (kb) 
^ l x (kb) K a (ka) J 



(11) 



For a gradient coil of length L, the Fourier expansion around the geometric center for the current of the inner coil in a 
self -shielded design is: 

(z) =£ j n a sink n z for|z|*-| (12), 

n=l £ 



j^(z) = 0forlzl>| (13), 

where j n a are the Fourier coefficients of the expansion and sin k n z represents the antisymmetry condition of the current 
around the origin. For a coil of length L, the current is preferably restricted to become zero at the ends of the coil which 
restricts the values that k n can take. Thus, the allowable values for k n are: 

a 



\; ( Z ) = 0forlzl = |^k n = ?!p 



The Fourier transform of j<f, a (z) is: 



j«r(m,k> =£iij n a n/ n (k) 

n=l ^ 



(15), 



with 



4r (k) = 



-sin(k-k n )| sin(k + k n )| 



Due to the symmetry requirements, the dependence of ^(k) in the variable k is: 

f n (-k)-f„ (k) 



(16) 



(14). 



(17). 
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From the expression of the Fourier component of the current, the expressions for the gradient magnetic field B z 
the stored magnetic energy W in the gradient coil, and the z-directed thrust force F z are: 



10 



15 



B 



2n 



f-dlcsirikz k£ »j n >„(k) 

J ~°° , n=l 



I 0 (kp) K t (ka) 



i _ l 1 (ka) Kj (Kb) 
l x (kb) K x (ka) j 



(18), 



W = - 



f~dkl 1 (ka) K 1 (ka) 



^_ I x (ka) K x (kb) 
^ I x (kb) K x (ka) J 



E E j n a jn a . (k)Hr n . (k) (19) . 



n-1 n'=l 



20 



25 



F z net =aE i^j n a /-dk;|; n (k) 

n=l ^ J — 



(ka) . 
B po< k > -TS^7 B P0(k) 



(20) . 



I x (kb) 

Because B p0 a ' b (z) is an odd-symmetric function with respect to z, its Fourier transform is defined as: 

Bp 0 ' b (k) = -i f yz (sinkz) B p a 0 ' b (z) 

= -iB p a ' b (k) (21a) , 



30 with: 



B p a 0 ' b (k) =[yz (sinkz) B p a ' b (z) 



(21b) 



55 Equation (20) then has the form: 

L 



40 



45 



50 
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(22) 



From the expressions for the magnetic field, the stored energy, and the z-directed thrust, the functional E can be 
defined as: 

E(j„ a )=W-f: ^(B z ( ?j )-B zsc (f ))- N £ ^(f^-F^) (23), 

j=l J j-N^l V ' 

where B zSC and F zSC net are the prespecified constraint values of the magnetic field at the constraint points and the 
z-directed thrust force, respectively. 

Minimizing E with respect to j n a , one obtains a matrix equation for the j n i a as follows: 



i _ I 1 (ka)K 1 (kb) 



E j^|-^/_"dkI 1 (ka)K 1 (ka; 

= ~22\ dkksinkz.I n (kp.) K. (ka) 1 - — 1 1 

JJ - 1 0 1 1 I. ( kb) Kj (ka) ) 



i|L(k)* n . (k) 



i^kbJK-tka) r n v 



(24a) , 
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10 



or 

N. +1 



E jn"C n , n = E A D (24b) 

n'=l j=l J J 



Depending on the value of j, the D jn matrix can correspond to either the magnetic field expression or to the z-directed 
thrust force expression. 

Specifically, the conversion for the D jn matrix is: 



^jn | r net 



B-field for j =1,*^ (25) 
F 2 net for j =N 1 +1, N x +1 



15 For the expression of B p0 a ' b (z), data is obtained from the actual main magnet 10. Due to the antisymmetry of the 

B p0 a ' b (z), only positive z-values are generated. In order to incorporate the generated data for the radial component of 
the magnetic field into Equation (24a), consider the following: 

B p a o b (z) =E B Y a ' b (z )5(z-z y ) (26), 

20 Y =0 Y Y 

where Q corresponds to the total number of points along z, and B y a - b is the corresponding value of B p0 a b (z) at the 
location z y . The Fourier transform B p0 a ' b (k), is: 

B p a 0 ' b (k) =2E B y a ' b (z Y )sinkz Az Y (27). 

25 P Y =0 Y Y Y 

Furthermore, the generated values for B p0 ab (z) are for the axial distances which are larger than the half length of 
the gradient coil and up to the point where there is no significant action between the current density of the gradient coils 
and the main magnet's radial components. Truncating the infinite summations at M terms, the matrix representation for 
Equation (24b) becomes: 

30 

J a C=A. D or J a = X D C (28) 

where J a is a 1xM matrix, C is a MxM matrix, X is a IxN-, +1 matrix, and D is an N-j+lxM matrix. 

Using the constraint equation for the gradient magnetic field B z and the z-directed thrust force, the Lagrangian 
multipliers are determined. The Fourier components of the current density for the primary or inner gradient coil in matrix 
35 form are: 

/ = § z [2 C"V]" 1 D C" 1 (29). 

Finding the expression for the Fourier components of the current for the inner coil and with the help of Equation 
(14), the continuous distribution of current density of the inner coil is generated. To determine the current for the outer 

40 coil, the shielding relationship between the Fourier transform for both current densities of Equation (8) is used. An inverse 
Fourier transform is then used to obtain the continuous current distribution for the outer coil. 

The next step is the process for discretization of the continuous current distribution for both coils. The continuous 
current distribution is divided into positive and negative current regions. Integrating the area underneath each region, 
the total current contained in each region is obtained. When the current for all of the regions of the cylinder is calculated, 

45 discrete current loops are positioned on a dielectric former to mimic the behavior of the continuous current pattern. Each 
region is filled with discrete wires carrying the prescribed amount of current. The amount of current is the same for each 
wire loop, in the preferred embodiment. In each region, the continuous current density is divided into smaller segments 
which correspond to the selected equal amount of current. The selected current amount may be iteratively adjusted in 
order to match the continuous current densities in the selected current regions more closely. Each wire is placed at the 

50 midpoint of the corresponding segment in order to obtain an equal distribution from both sides of the segment. This 
current distribution is then analyzed to calculate the generated magnetic field to assure that the desired magnetic field 
is, in fact, obtained. Alternately, one can use the center of mass scheme taught in previously referenced U.S. Patent 
No. 5,296,810. 

With continuing reference to FIGURE 2, in the preferred embodiment, the inner or primary coil has a radius of about 
55 340 millimetres and a length of about 700 to 900 millimetres. The outside or shield coil is about 385 millimetres in 
diameter. For the Fourier series expansion, ten points provide a reasonable definition of the current density, although 
larger or smaller numbers may be chosen. Five constraint points are preferably chosen to define the characteristics of 
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the field and the thrust force. The field is to be constrained inside a 40 cm. to 50 cm. generally ellipsoidal imaging volume. 
The first constraint point establishes a gradient field strength, e.g., about 1 3.5 mT/m. The second constraint point defines 
the linearity of the field along the gradient axis. Preferably, at the edges of the 25 cm. dimension of the volume, the 
magnetic field is confined to vary not more than 5%. The remaining two constraints define the uniformity of the field 
across a plane perpendicular to the gradient axis. Preferably, the magnetic field is constrained to stay within 7% of its 
actual value at a radial distance of about 22 cm. from the center of the coils. 

The last constraint point defines the value of the z-directed thrust force which is preferably less than -1 .Oe -08 New- 
tons. A suitable primary gradient coil meeting these constraints is illustrated in FIGURE 2A and a suitable outer or shield 
gradient coil meeting these conditions is illustrated in FIGURE 2B. It will be noted that the primary gradient coil has 
several force correcting windings 64 of reversed polarity at its ends or extremes. If these reverse polarity windings were 
removed, the gradient coil would suffer a net longitudinal force component in excess of the above-discussed constraints. 

Although the above-described method calculates an ideal current distribution, the force correcting current for can- 
celling the axial force can be determined iteratively. More specifically, several coil windings are disposed at the ends of 
one of the primary and secondary coils, preferably the primary coils. The net force on the gradient coil assembly during 
a gradient current pulse is measured. The current flow through the force correcting windings (or the number of force 
correcting loops) is iteratively adjusted until the axial force is substantially zeroed. 

With reference to FIGURES 3Aand 3B, the gradient coil assembly further includes an x-gradient coil and ay-gradient 
coil. The x and y-gradient coils are of substantially the same construction, but rotated 90°. Of course, because one is 
laminated over the other, they will have a small difference in radius. The x and y-gradient coils each include four sub- 
stantially identical windings arranged symmetrically on either side of the isocenter. Each of the quadrants of the primary 
x or y-gradient coil contains a winding substantially as illustrated in FIGURE 3A. Each of the quadrants of the shield x 
or y-gradient coil contains a winding substantially as illustrated in FIGURE 3B. A current distribution is again calculated 
to produce the selected gradient strength, minimize stored energy, and zero the net lateral Lorentz forces. 

For the transverse x or y-gradient coils, the total current density can be represented as: 



where 5(p-p 0 ) is the restriction that the current is confined on a cylindrical surface of radius p 0 . Again, the x and y-gradient 
coils are self-shielded. That is, they have an inner or primary coil of radius a and an outer secondary or shield coil of 
radius b. The current density for the outer coil of radius b is analogous to Equation (30). The interaction between the 
components of the current density and the corresponding components of the main magnetic field yields an x-directed 
Lorentz force for self-shielded x-gradient coils and a y-direction force for self-shielded y-gradient coils. The following 
discussion focuses on x-gradient coils. However, it is to be appreciated that the same discussion is equally applicable 
to y-gradient coils which will be of the same construction but rotated 90° and of slightly different radius. 

The Lorentz equation describes the elemental force experienced by a current element Id1 in the presence of the 
magnetic field B 0 : 



J a (f) = j z)<V j z a (0,z)(3 z 5(p-p 0 ) 



(30) , 



dF = Idl x § 0 



Ola) , 



or 



dF = (j*'" ( z) cos (<|>) <X 0 + j a ' b { z ) sin «t») <3 2 )dz (a, b) d<|> 




(31b) , 



where B p0 ab (z) and B z0 ab (z) are radial and axial components of the main magnetic field B 0 at radii locations (a,b). 
Integrating Equation (31b) results in a net force along the x-axis: 



F" et = F x n 



net 




(32a), 



with: 



F x net = 4naJ o ~{j 0 a (z)B z a o (z) - j z a ( z ) B p a 0 ( z ) > dz 

+ 4nb/ o ~{j^ b (z)B z b (z) - j z b (z)B p b 0 (z)}dz (32b), 



where B p0 a ' b (z) and B z0 a ' b (z)are known real functions of z. The Fourier transform pairs for j^' 6 , j z a ' b , B p0 
are: 



a ' b , and B z0 a ' b 
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j; ,b (k) = /_;V b (z)e- ikz dz 
j, a ' b (z) = ^L|_-j^ b (k)e- k Mk 
j a ' b (k) = /~j a ' b (z) e - ikz dz 

1 f +«,_. a,b , _ +ikz^^. 



j z a ' b (z) = 



2n 



/~j a ' b (k)e 



B a o' b (k) = /jB p a 0 ' b (z)e- ikz dz 



(33a) , 
(33b) , 

(33c) , 
(33d) , 

(33e) , 



1 f "K 0 ' b (k)e tikz dk 



(33f), 



B z0 ' b (k) = jj-B z a 0 ' b (z)e- ik Mz (33g) , 



B a 0 ' b (z) = ^- f "B z a o b (k) e +ikz dk (33h) , 

2n J -~ 



Using the constraint: 



V • J = 0 



and the shielding conditions, the following relationships are derived: 

. b a l 'l( ka ) . a 

The expression for F x net from Equation (32b) is rewritten as: 



F x net =2a/ o ~j 0 a ' b (k)| 



I\(ka) 



I'^kb) 



B z *o (k) ^-1-B; a (k) 



B z * b (k) + -L Bp * b (k) 



1_ 

kb" 



(35), 



(34a), 
(34b), 



where * represents the complex conjugate of the Fourier transforms of the two components of the main magnetic field. 

For a real magnet design, there is no useful analytic expression forthe B z0 a ' b (z) and B p0 a b (z). In orderto incorporate 
analytically the numerically generated data for both components of the magnetic field along the central z-axis of the 
magnet, consider the following representations: 
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B 



a,b 
zO 



(k) 



-2± 

Y=0 



B z a ; b (z Y ) cos kZyAZy, = B z * 0 (a ' b) (k) 



(36a) , 



B 



a,b 
PO 



(k) 



>.± 



= 22. B^Nz^sinkz^-B^ 3 '^ (k) 



*(a,b) 



(36b) , 



where B z0 a ' b (z) is symmetric around z and B p0 a ' b (z) is antisymmetric around z. 

The restriction to the inner coil length, the confinement of the current density to a cylindrical surface, the azimuthal 
and axial symmetries for the and j z a and the demand that the current density obeys the continuity equation, provides 
a Fourier series expansion for both components around the geometric center of the coil: 



j 0 a ((t>,z) =cos<{>X) j^cosk n z for |z| ^ 



n=l 



(37a), 



j z (4>,z) =sin 



*£4 



n-i k n a 



sink n z for |z| * 



(37b) , 



where j^ n a are the Fourier coefficients, L represents the total length of the inner coil, and k n =2nn/L because the current 
cannot flow off of the ends of the cylinder. Furthermore, both current components are zero for lzl>L/2. 



B 



2n m ' 



E e* im *f-dkke ik % a (m,k) 



l m (kp)K' m (ka) 



1 - 



I' m (ka)K' m (kb) 
I' (kb)K' (ka) 



(38), 



The general expression for the magnetic field for a self-shielded gradient coil in terms of the Fourier transform of 
the current density is: where j (() a (m, k) is the double Fourier transform of a (c|),z). Since the azimuthal dependence of 
\q a is proportional to cos(c|)), the Fourier transform of \^ a is non-zero when m=±1 . In this case the two-dimensional Fourier 
transform of the current density is: 



j<T(±i,k) =-iE^j<,X(k> 

Z n=l Z 



with: 



(39), 



sin(k-kj| sin(k + k n )| 



(40), 



where >F n (k) is an even function of k and a (+1 , k)=j ( | ) a (-1 , k). Therefore, the expression of the gradient field has the form: 

Po aL 



B = 
z 



4n 

I, (kpJK'^ka) 



cos<|)E j^/jikcoskzt^ (k)k 
^ I'jtkaJK'^kb) N 



^ I'^kbJK'^ka) J 

In a similar fashion, the stored magnetic energy in the system is: 



(41) 
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a. L ^ °° 00 
w = - ° E E j*nj/n-/"dkUf n (k)Hf n , (k) 

1 ° n=l n* =1, J -~ 

I'^kaJK'^kb) x 



I\ (kaJK'^ka) 



1 - 



I'^kbJK'^ka) J 

and the expression for the x-directed Lorentz force is written as: 



(42) , 



F x net = aL£ j^/~t|; n (k)dk 

n=l J 0 

I'itka) 



l\(kb) 



B z *o(k) + -^-B p * 0 a (k) 



B z * 0 b (k) + -±-B p * 0 "(k) 



1 p »b 

kb 



(43) . 



The functional E is constructed and has the identical form of Equation (22). However, in the present case, the 
dependence of E is over j^ 8 and the force is directed along the x-axis instead of the z-axis. Thus, minimizing E with 
respect to j^ 3 yields a matrix equation for a as follows: 



Hf n (k)tj/ n , (k) 



v> . a J aLn /•- I 'i (ka)K\ (kb) 

E j/n-i^f dkl' (ka)K' (ka) 1- — ! — 

n-=i * [ 2 io 1 1 ^ I'^kbJK'^ka) j 

^ „ ( I', (ka)K'. (kb) 

= 2^cos (<J>.) / dkkcoskz.I.. (kp. ) K', (ka) 1- — 1 3 

j-i 3 3 J ° 3 1 3 1 { I'^kbJK'^ka) t 

"£ x,25/-dkt n (k){ 



B z0 (k) + 



(ka) 



B p0 (k) 



B zo< k > + ^B p * 0 b (k) 



kb 



M N^l 

E j-t>n' C l>>n = E ^jD jn 
n'=l j=l 



(45), 



I'itkb) 

Truncating the infinite summations at M terms, and using compact notation, Equation (38) simplifies to 
or: 

J a C=A D or J a =A D C" 1 

where J a is a 1 xM matrix, C is a MxM matrix, X is a 1 xN-, +1 matrix, and D is a +1 xM matrix with 

(47) . 

for j = N x + 1, N x + 1 



(44) 



(46), 



J B- field for j = 1,^ 



Because Equation (46) is the same as Equation (28), the expressions of the continuous current distribution for both 
coils is found by following the steps described above. To discretize both current densities, one first considers the con- 
tinuity equation for the current density: 

V * J = 0 (48) . 

In analogy with the magnetic field, where a vector potential is introduced, the current density is expressed as a curl of 

the stream function s as: 

J=V*S (49). 
Because the current is restricted to flow on the surface of a cylinder of radius a=p a and has only angular and axial 
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dependence, the relation between the current density and the stream function in cylindrical coordinates is: 

a a a a ^Sn A 1 A 

j + (♦,z)o^+J 2 (*,z)a z = -^a ri -^a z (50), 

and S p is found from: 

S p «t>,z) = -aj*_*d<t>' j z a (<D',z) (51). 
The contour plots of the current density are determined by: 

S p (x,z) =( n-lj S^+S^ for n = 1, N (52), 



15 the amount of the current between two contour lines. The determination of S inc is 



where N is the number of the current contours, S mjn is the minimum value of the current density, and S jnc represents 

mil 

S - S ■ 

r> max mm /cn\ 
S inc = M 

cun 

20 follow the flow of the current and the distance between them corresponds to a current equal to an amount of S inc in 
amps. The discrete wires are positioned to coincide with these contour lines. 

In the preferred embodiment, the self-shielded force free x-gradient coil has a radius of the inner coil of about 1/3 
metre and a length which is about -100 to 900 mm. The radius of the outer coil is about 0.4 metres. For the design of 
this self-shielded coil, the number of Fourier coefficients which are used is equal to the total number of constraint points 

25 increased by 1 . In addition, four constraint points are chosen to specify the quality of the magnetic field inside a 40-50 
cm. by 25 cm. ellipsoidal volume, and to eliminate the net Lorentz force for the whole gradient coil system. The first 
constraint sets the strength of the gradient field, e.g. 13.5 mT/m. The second constraint specifies the linearity of the 
gradient field along the x-axis and up to the distance of 25 cm. from the isocenter of the gradient field, e.g., 5%. The 
third constraint specifies the uniformity of the gradient field along the z-axis at a distance up to +6 cm. from the isocenter 

30 of the gradient field, e.g., 1 2%. The fourth constraint specifies the value of the net x-directed force on the gradient sets 
location, e.g., 1x10 -7 Newtons. 

The presence of this set of constraints generate a continuous current distribution for both the inner and outer coils 
which, when discretized, takes the form illustrated in FIGURES 3A and 3B. More specifically, applying these design 
requirements generates a continuous current distribution of j za and j zb . Only the z-component of the current density is 

35 necessary for the creation of a discrete current pattern. The constraints are also chosen such that an integer number 
of turns for the inner and outer coils with a constant amount of current per loop can be obtained. Using the stream 
function technique, the discrete current distribution of the inner coil for the lateral force free configuration is obtained as 
illustrated in FIGURE 3A. It will be noted that the inner or primary coil has a thumbprint winding 70 and an additional or 
force offsetting winding 72. The current in the force correcting winding effectively flows in an opposite spiral direction to 

40 the primary coil portion 70. The current distribution and number of turns in the force cancelling coil portion 72 is selected 
to zero the net x or y-direction force generated by the primary coil section 70 of FIGURE 3A. Depending on the nature 
of the B0 field components in a particular design, similar force free results can also be obtained by placing a winding 
having an opposite sense than the secondary x or y-gradient coil near an outer edge of that coil. The amperage through 
the force offsetting windings or the number of windings are adjusted until a current density is achieved which produces 

45 an x or y-force component which is equal and opposite to the force component generated by the other x or y-gradient coils. 

Although force balancing coils are applied to only the primary winding of shielded gradient coils in the preferred 
embodiment, it is to be appreciated that the primary and secondary coils can each be individually force balanced. Further, 
it is to be appreciated that the thrust force balancing is also applicable to gradient coils which are not self -shielded, i.e., 
which have only a primary winding and no secondary shield winding. 

50 



Claims 



1 . A magnetic resonance imaging apparatus which includes a main magnet (1 0) for generating a temporally constant 
55 magnetic field through a patient receiving bore (12) and a gradient field coil assembly (30) and an RF coil (40) 

around the bore (12) and in which the temporally constant magnetic field has components which interact with current 
pulses applied to principal windings (32, 36) of the gradient field coil assembly (30) to generate magnetic field 
gradients in such a manner that a net force is experienced by the gradient field coil assembly (30), characterised 
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by force correction windings (64, 72) mounted adjacent edge areas of the gradient field coil assembly (30) which 
carry currents in an opposite sense to the principal windings (32, 36) such that the net force is substantially cancelled. 

A magnetic resonance imaging apparatus as set forth in claim 1 wherein the patient receiving bore (1 2) has a length 
to diameter ratio of less than 1 .5:1 . 

A magnetic resonance imaging apparatus as set forth in claim 1 wherein the patient receiving bore (1 2) has a length 
to diameter ratio of substantially 1:1. 

A magnetic resonance imaging apparatus as set forth in any preceding claim wherein the gradient coil assembly 
(30) is a self-shielded gradient coil having a primary coil (32) and a shield coil (36), current pulses passing through 
the primary and shield coils (32, 36) generating magnetic fields that combine within the gradient coil assembly (30) 
to form a linear magnetic field gradient and tending to cancel outside the gradient coil assembly (30), and the current 
pulses that pass through the shield coil (36) interact with radial magnetic field components to generate a first force, 
the primary gradient coil (32) being wound such that the current pulses passing through the primary gradient coil 
(32) interact with the radial magnetic field components to generate a second force that is substantially equal and 
opposite to the first force. 

A magnetic resonance imaging apparatus as set forth in any one of Claims 1 to 3 wherein the gradient coil assembly 
(30) includes primary and secondary gradient coils (32, 36) which each include four symmetrically arranged thumb- 
print coils (70) for generating a magnetic field gradient transverse to a longitudinal axis of the bore (1 2), the thumb- 
print coils (70) further generating a net transverse force, the force correction windings including force offsetting 
current loops (72) disposed adjacent an end of at least one of the primary and secondary gradient coils (32, 36) for 
carrying current flows which circulate in an opposite direction to a most adjacent thumbprint coil (70) for offsetting 
the net transverse force along the transverse direction. 

A magnetic resonance imaging apparatus as set forth in any one of Claims 1 to 3 wherein the gradient coil assembly 
(30) includes a primary z-gradient coil (32) including a series of loop coils (60) which are mounted around a cylindrical 
former (34) for generating a z-gradient field and the force correction windings include loops (64) extending around 
ends of the cylindrical former (34) and connected in series with the z-gradient loop coils (60) for carrying common 
current pulses but in an opposite direction. 

A magnetic resonance imaging apparatus as set forth in any one of Claims 1 to 3 wherein the gradient coil assembly 
(30) includes a primary z-gradient coil (32) and a secondary, shield z-gradient coil (36) and the secondary z-gradient 
coil (36) is disposed surrounding the primary z-gradient coil (32), current pulses applied to the primary and secondary 
z-gradient coils (32, 36) causing magnetic fields which combine within the gradient coil assembly (30) to generate 
a linear gradient axially along the coil assembly (30) and which generally cancel exterior to the coil assembly (30). 

A magnetic resonance imaging apparatus as set forth in claim 7 wherein the primary and secondary gradient coils 
(32, 36) include annular loops (60, 62) and the force correction windings include loops (64) disposed on the ends 
of at least one of the primary and secondary coils (32, 36). 

A magnetic resonance imaging apparatus as set forth in any one of Claims 1 to 3 wherein the gradient coil assembly 
(30) has a z-gradient coil including: a cylindrical former (34); a first plurality of loops (60) extending around the former 
(34) on a first side of the isocenter of the temporally constant magnetic field, the first plurality of loops (60) carrying 
current pulses primarily in a first direction circumferentially around the former (34); a second plurality of loops (60) 
extending around the former (34)on a second side of the isocenter mirror image symmetric about the isocenter with 
the first plurality of loops (60), the second plurality of loops (60) carrying the current pulses primarily in a second 
direction circumferentially around the former (34), the second direction being opposite to the first direction, the 
current pulses passing through the first and second plurality of loops (60) interacting with the radial components of 
the temporally constant magnetic field to cause a net axial force; first force offset coils (64) extending around the 
former (34)on the first side of isocenter displaced from the isocenter, the first force offset coils carrying the current 
pulses in the second direction; and second force offset coils (64) extending around the former (34) on the second 
side of isocenter and displaced from the isocenter, the second force offset coils (64) carrying the current pulses in 
a first direction and the current pulses passing through the first and second force offset coils (64) interacting with 
the radial components of the temporally constant magnetic field to cause an axial force which offsets the net axial 
force. 
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10. A method of designing gradient coils for a magnetic resonance imaging apparatus in which coil windings patterns 
are constrained to have a preselected gradient field strength, at least a preselected gradient field linearity and, at 
least a preselected gradient field uniformity, characterized by: constraining the coil windings such that the gradient 
coil assembly experiences substantially no axial or transverse force when current pulses flowing through the gradient 
coil windings interact with radial and axial components of a main magnetic field within which the gradient coil is 
disposed. 
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